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Four-loop Standard Model effective potential at leading order in QCD 
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The leading QCD part of the four-loop contribution to the effective potential for the 
Standard Model Higgs field is found. As a byproduct, I also find the corresponding 
contribution to the four-loop beta function of the Higgs self-interaction coupling. 


Contents 


I. Introduction 1 

II. Three-loop and four-loop integral basis 3 

III. Effective potential in terms of bare quantities 5 

IV. Effective potential in terms of renormalized quantities 6 

V. Discussion 11 

References 12 


I. INTRODUCTION 

The effective potential [1-3] is an important tool for analyzing spontaneous symmetry breaking 
associated with scalar field vacuum expectation values (VEVs). In the Standard Model, it provides 
a quantitative link between the Lagrangian parameters and the VEV of the Higgs field. The 
fact that the Higgs boson mass is near 125 GeV implies that the electroweak vacuum is close to 
metastable, motivating a program of precise study of the stability criteria [3-21]. Of more general 
importance is the fact that the effective potential minimization condition allows one to determine 
and eliminate one of the Lagrangian parameters of the theory, typically the negative Higgs squared 
mass parameter, in favor of the radiatively corrected VEV. 

The effective potential V e s((f) can be obtained as the sum of one-particle irreducible vacuum 
Feynman graphs, computed in terms of particle masses and couplings that depend on a constant 
background scalar field (f). In the normalization conventions of the present paper, the canonically 
normalized Standard Model Higgs complex doublet field has a tree-level potential 


V = m 2 <h t <h + A(<h t ^>) 2 , 


( 1 . 1 ) 
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where A is the Higgs self-interaction coupling, and the negative Higgs squared mass parameter is 
m 2 . The real neutral part of $ is given by (cj>-\-h)/\/ 2, where 4> is the constant background field and 
h is the physical Higgs real scalar boson field. The complete set of 1-loop and 2-loop contributions 
to the effective potential in Landau gauge are known for the Standard Model [22] and for a general 
renormalizable field theory [23]. Also known [24] are the 3-loop contributions that only involve the 
strong coupling g 3 and the top-quark Yukawa coupling y t . Contributions from Goldstone bosons 
can be resummed [25, 26] in order to avoid potential infrared singularities and spurious imaginary 
parts. The value of the background field at the minimum of the effective potential is the radiatively 
corrected VEV of the Higgs field. 

The purpose of this paper is to extend the existing calculations of the effective potential V e s (</>) 
by obtaining the 4-loop contributions that are leading in the strong coupling 53 , using dimensional 
regularization [27-31] and the MS renormalization scheme [32, 33]. These contributions come 
from those diagrams that involve only quarks, gluons, and QCD ghost fields. I will work in the 
approximation that all quarks are massless except the top quark. This is an excellent approximation 
beyond 1-loop order, due to the small magnitudes of the Yukawa couplings of the bottom and 
other quarks. Then, in dimensional regularization, at least one top-quark loop must be present in 
a diagram in order for the contribution not to vanish. At loop order £, the resulting leading QCD 
contribution is proportional to ^ t 2 multiplied by a polynomial of order £ in ln(f), where 

t = yffi/2, ( 1 . 2 ) 

is the field-dependent MS top-quark squared mass, and 

In (f) = In (t/Q 2 ), (1.3) 

where Q is the MS renormalization scale. 

The organization of the remainder of this paper is as follows. In section II, I review the basis 
of scalar integrals used in the calculation. The effective potential in 

d = 4 — 2e, (1.4) 

spacetime dimensions is given in section III in terms of bare quantities and the basis integrals. In 
section IV, the bare parameters are re-expressed in terms of MS quantities to obtain the effective 
potential in that renormalization scheme, after expanding in e. (This is more efficient than doing 
a separate calculation of counterterm diagrams.) In the process, I obtain the leading QCD con¬ 
tribution to the 4-loop beta function for A, from the requirement that poles in e do not appear in 
the effective potential when written in terms of the renormalized parameters. Some of the results, 
when given in general form in terms of group theory invariants, are rather lengthy and therefore 
are provided in ancillary electronic files rather than in print. Section V concludes with some brief 
comments on the numerical impact of the new results. 
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II. THREE-LOOP AND FOUR-LOOP INTEGRAL BASIS 

In the approximation of this paper, the only mass scale (other than the renormalization scale) is 
the top-quark mass. Therefore, it is convenient to write results in terms of Euclidean momentum 
integrals with each propagator having dimensionless mass 0 or 1. The dependence on the bare 
top-quark mass is then restored by dimensional analysis. Only integrals having an even number of 
massive propagators meeting at each vertex are needed in this paper. Momentum integrations in 
d dimensions are normalized by 



so that the 1 -loop vacuum master integral is defined by 


( 2 . 1 ) 


A = 


M T 


T(1 — d/2) 
( 47 r ) rf / 2 


( 2 . 2 ) 


At 2-loop order, no new master integral appears. The necessary 3-loop and 4-loop integrals have 
been studied and used in refs. [34-51]. Important applications include the calculations of the 4-loop 
QCD corrections [41, 45, 46] to the p parameter and decoupling rules for as and light quark masses 
across heavy quark thresholds [43, 51]. Figure 2.1 shows a basis for the master integrals needed [39] 
for single-scale gauge theories at 3-loop order and 4-loop order. Each solid line represents a massive 
propagator denominator, and each dashed line represents a massless propagator denominator, and 
the Euclidean loop integrations are normalized according to eq. (2.1). So, for example, 


142 J P JqJkP 2 q 2 {k 2 + l)[(p + q + k) 2 + l]' ^' 3) 

Also needed in the basis are products A 2 , A 3 , A 4 , A/ 40 , and A/ 42 . All of the integrals used in this 
paper are reduced to the basis by repeated application of the integration by parts method [52], 
using a strategy similar to that described in ref. [39]. 

The integrals 1. 42 , I 53 , and have non-zero masses confined to a single 1-loop self-energy 
subdiagram, and are therefore known analytically in terms of T functions. In general, it is sufficient 
to have results for the basis integrals as expansions in e. However, with the basis chosen here, the 
coefficients of the basis integrals have poles in 1 /e in addition to the poles inherent in the basis 
integrals.^ This means that it is necessary to have the expansions to certain positive powers 
of e in most of the cases. The coefficients of expansions in e of the other integrals have been 
given numerically with high precision and to sufficiently high order in e in [40], using the Laporta 
difference equation method [53]. In principle this is enough for practical purposes, but it is nice 
to have analytical versions as well. These have been provided in refs. [40, 47, 49, 50]. Table 2.1 
shows the order in e to which each basis integral is needed in the present paper, as well as the 
highest order to which it is known analytically in terms of simple e-independent sums, and the 


t For an alternative basis with the nice property that coefficients do not contain extra poles in e, see ref. [44, 48]. 
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FIG. 2.1: The 3-loop and 4-loop scalar basis integrals with one mass scale and an even number of massive 
propagators at each vertex. Massive propagator lines are solid, and massless propagator lines are dashed. 
The first subscript is the total number of propagators, and the second is the number of massless propagators. 


TABLE 2.1: Summary of present analytical knowledge of 3-loop and 4-loop basis integrals needed in this 
paper and depicted in Figure 2.1. The second row shows the order in the expansion in e needed here. The 
third row shows the highest order in the e expansion to which the integral is presently known analytically 
in terms of simple e-independent sums, while the fourth row gives the source reference for that expansion. 
All integrals were also previously evaluated numerically to the necessary orders and beyond in ref. [40]. 


integral 

ho 

ho 

hi 

hi 

ho 

h2 

h3 

hi 

I 72 

ho 

hi 

ho 

needed 

e 3 

6 4 

e 4 

e 3 

e 2 

e 3 

e 2 

e 1 

e 1 

e 1 

e° 

h 

known 

e 5 

e°° 

e 4 

6 °° 

e 2 

e 3 

e 5 

e°° 

e 3 

e 1 

e° 

e 1 

source 

[49] 

[40] 

[50] 

[40] 

[50] 

[49] 

[47] 

[40] 

[49] 

[50] 

[50] 

[50] 


source reference that provides that expansion. The (probably) transcendental numbers appearing 
in these coefficients up to the orders needed in this paper are tv, ln(2), and 


Cn 





k=l 


Li n (l/2) = ^ 
k=1 w 


oo n 


EE 


n 5 k ’ 


(2.4) 

(2.5) 

( 2 . 6 ) 


although the last quantity cancels out of the results below. (The absence of this quantity could 
presumably have been made manifest by using the alternative basis of [44, 48].) 

































5 


III. EFFECTIVE POTENTIAL IN TERMS OF BARE QUANTITIES 

In this section, I find the 4-loop effective potential in terms of the bare quantities in d = 4 — 2e 
dimensions. These include the external scalar field (f>B and the bare Yukawa coupling yts and 
QCD coupling g^s- In the next section, the results will be converted to MS parameters. The loop 
expansion for the effective potential is written as 


OO 

r=» = E<> 


e=o 


(3.1) 


The tree-level potential is 



m 


B ,2 


+ T^» 


(3.2) 


where A b and m 2 B are the bare Higgs self-coupling and squared mass parameter, respectively. The 
latter will play no role in the following. 

At each loop order, the contribution to the effective potential is given by the sum of 1-particle 
irreducible Feynman diagrams with no external legs and containing only quarks, gluons, and QCD 
ghosts. The pertinent contributions at loop order £ > 1 are proportional to 4 ^ 2 i where 


tB — Vw^b/Z 


(3.3) 


is the bare field-dependent top-quark mass. Results below will be given in terms of group theory in¬ 
variants: the dimension of the fundamental representation N c , the Casimir invariants of the adjoint 
and fundamental representations Cq and Cf, the Dynkin index of the fundamental representation 
Tp, and the number of quark flavors n q . In the Standard Model, these are given by 

C G = N c = 3, (3.4) 

N 2 - 1 

C f = ^-=4/3, (3.5) 

T f = 1/2, (3.6) 

n q = 6, (3.7) 

but leaving them general provides more information for comparisons and checks. Diagrams at 
2-loop order and higher are calculated with a gluon propagator 

-i[g^/ P 2 -(i -£)pV7(p 2 ) 2 ], ( 3 - 8 ) 


where £ = 1 for Feynman gauge and £ = 0 for Landau gauge. The dependence on the (bare) QCD 
gauge-fixing parameter £ cancels at the level of the basis integrals, providing a stringent check. 
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The contributions involving only quarks, gluons, and QCD ghosts, up to 3-loop order, are [24]: 


V W 

V B 

y( 2 ) 

V B 

y(3) 

V B 


—AN c t d g 2 A/d, 

N c C F g 2 B t d B 2 A 2 {d - 1 )(d - 2)/(d - 3), 


gi B t 3d,2 ~ 4 N c C F {c G 


(2 ~df 


L2(d — 4) 2 (d — 3) 


+ 


+ 


+ 


+ 


I42 


+ C F 


(d — 2) 2 (2d 2 - 17d + 32) 

2(d — 4) (2d — 7) 

(d - 2) 2 (-d 5 + 13d 4 - 67 d 3 + 181d 2 - 274d + 188) 
2(d — 4) 2 (d — 3) 2 
(2 - d)(2d 3 - 21 d 2 + 67d - 68) 


3 (3-d) (d 3 - 13d 2 + 50d - 48) 

+ 4(d-4) 2 ' 

(d — 6 )(d — 3)(d 2 — 7d + 8) 

2(d — 4) 2 


*40 


A 3 


I42 


+ T f 


2(5 - d){d-2f 


(d — 4)(d — 3) 
d 3 — 7d 2 + 6d + 16 . 4(d — 3)(d — 2) 

(d — 6)(4 — d) Im + ( " a _ 11 -- h2 


7-2d 


L(d — 6)(d — 4)(d — 3) 


^ 3 


*40 


(3.9) 

(3.10) 


(3.11) 


For the 4-loop order contributions involving quarks, gluons, and QCD ghosts, there are 51 
Feynman diagrams, which are reduced to linear combinations of the 13 integrals from the set 


T — {A 4 , Af 40 , AI 42 , hi, I 53 , ho, I 02 , 


63) *64) 


I72, I73, hi, ^ 93 }) 


(3.12) 


using integration by parts identities. The four-loop effective potential contribution is then organized 
in terms of the group theory invariants from the set 

Q = { Cq , C G T F , CcT F n q , C G C F , Cp, C F T Fl C F T F n q , T B , T B n q , T B n 2 }, (3.13) 

so that the result is written as: 

Vb ] = 9l B t 2 B d - 6 C F N c J2Yl (3.14) 

Q X 

The 130 coefficients V b \q,T) are rational functions of the spacetime dimension d. Although 58 
of them vanish, this list of coefficients is still rather lengthy, so they are not shown in print here. 
Instead, they are provided in an ancillary electronic file called V4bare . txt included with the arXiv 
submission for this article. 


IV. EFFECTIVE POTENTIAL IN TERMS OF RENORMALIZED QUANTITIES 

In this section, I obtain the effective potential in the MS renormalization scheme by translating 
the bare quantities into MS quantities. Because f d d xV must be dimensionless in order to be 
exponentiated in the path integral, one must introduce an arbitrary regularization scale //, which 
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is related to the MS renormalization scale Q by [32, 33]: 


Q 2 = 47re lE p 2 . 


Then, in the MS scheme, one writes: 


4>b 

Z(f> 


XkB 


oo t 

‘+EE 


%n 


(l67r 2 Te n ’ 

=1 n =1 v ; 


n px k e (x k +^ 2 Y^ 


°e,n 


—' Z-j ( 16 iT 2 Ye n J 

= 1 n= 1 V ' 


(4.1) 


(4.2) 

(4.3) 

(4.4) 


The subscript B labels bare quantities, while the absence of a subscript B indicates the corre¬ 
sponding MS renormalized quantity. The exponent t is the loop order, while k is an index that 
runs over the list of Lagrangian parameters, including x k = AThe mass dimensions of the 
bare parameters determine that p\ = 2 and p g3 = p yt = 1 , in order that the renormalized couplings 
A, 53 , and yt are dimensionless and 4> has mass dimension 1. The counter-term quantities zf n and 
Zp k n are polynomials in the MS renormalized parameters Xj, and do not depend on e or <f). They 
are determined by the requirement that the full effective potential (and all physical observables) 
are free of ultraviolet poles in e when expressed in terms of the MS quantities. 

The anomalous dimension for 4> and the MS beta functions for the parameters x k are defined 
by 


7 

Pxk 


= -Q- 


dQ 


E=0 


dQ 


\ Q ik uz * ) ’ 


Q— 
* dQ 


n dx k 

E=0 ~ Q dQ + €pXkXk - 


(4.5) 

(4.6) 


Because the bare quantities 4>b and x k B do not depend on Q, the anomalous dimension and beta 
functions are determined by the simple pole counterterms, so that: 


7 


Px k 


OO 

V _-_7 W 

(167t 2 )^ 

OO ^ 

S (167T 2 /® ’ 


(4.7) 

(4.8) 


= ~f4,v 

tig = 


where the 7-loop contributions are: 


(4.9) 

(4.10) 



The higher pole counterterms are also fixed by consistency conditions 


£z tn 


t—n+1 


7 + 2 2^Px j Q ) z e-t',n-V 

3 J 


t—n+1 

E E<> 

i '=1 3 


dxj~ £ f/,n 1 


(4.11) 

(4.12) 


for l > n > 2. 

The coefficients zf and z* k n for i <3 are thus determined by the known results for the Standard 
Model beta functions and Higgs scalar anomalous dimension given in [54-60]. (Extensions to QCD 
4-loop and 5-loop order can be found in [61-66].) Keeping only the contributions needed for the 
approximation of the present paper, they are: 


* 1,1 = -N c yf + ..., (4.13) 

*2,i = 9ht(-2N c C F ) + ..., (4.14) 

* 2,2 = alvt^NcCp) + ..., (4.15) 

* 3,1 = 93 VtN c C F (^ 8(3 - Cq + -16C3^ Cf + ^16 + — n^j T F +•••, (4.16) 

* 3,2 = aiyt N cCF (24 Cq + 10CV — 16Tpnq/3) + ..., (4-17) 

* 3,3 = 9 ty$N c C F (-22C g /3 - 24 C F + 8T F n q /3) + ..., (4.18) 

* 1,1 = alyti-zCp) +..., ( 4 . 19 ) 

zf,! = 9 ivtC F (~C G - \c F + b -Tpn q ) + ..., (4.20) 

* 2,2 = aiytCp (y C G + \c F - 2 Tpn q ) + ..., (4.21) 

*h = alvtCp [~^g + f CgCf - T c r + (|f + 16 ^) CgTfn « 

+ (y - 16Cs) C F T F n q + ^ T 2 F rA + • • •, (4-22) 

* = CpC G + lcl-^C G Tpn q 

—^C F T F n q + —Tpn 2 J+..., (4.23) 

* 3,3 = gtytCp^-^fcl + y C G T F n q - y C G C F - \c 2 p + 6 C F T F n q - y^) + ..., (4.24) 

*T,i = 9 i(- l j-C G +‘rT F n ( i), (4.25) 

*fq = 93 ( — y C G + g CcTprig + C F T F n q ^j + ..., (4.26) 

* 2 % = 93 (yj“Q? — 3 - C G T F n q + -Tpnfj + ..., (4.27) 


while the zf n do not contribute at all at leading order in QCD. Now, expanding eq. (3.1) with 
eqs. (3.2), (3.9), (3.10), (3.11), and (3.14) to order 1/e, and requiring the 4-loop simple pole terms 
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to cancel, I find: 


z 4,i = yUi N cC F 


eg ( f C3 - 130C, C G T F + 886 - 40< 5 


^ m _ 661 16 t r 4 \ _ /896 826 . 1on . 4 t r 4 

+C G T F n, ( <3 - ^ + -45- J + CcC F (-3-g-Cs + 180 & + 77 

+CJ. ( 12Ca + 40( 5 till + lid j - 1SCfTf + C F T F n, - 111 


6 


64 o 2 /2728 32 A 

TpTlq ( -77777-7-C3 


6 


9 


V 


243 


+ 


z 4,2 = Vt93 N C C F 


f7m_44^ 88 _ 


LV 54 


3 J 


128Cs H38^ ^ 

-3 - W) CGTFn « 


+ ^-g- + ~3"^ 3 ) ^gC f + ( 2 + 48 Cs) Cp - 48CpT F 
451 112 .'i „ m . 32 m2 80 2 2 


+ ( —^ C 3 ) C F T F n q + —T F n q - 7^T F n q 


+ ..., 


z 4,3 


z 4,4 


yfglN c c F 
32 2 2 " 

~~Q~ T F n q 

ytatNcCp 

16 m9 2 1 
+ — T F n q 


2 322 562 2 146 

—61 C G + ——C G T F n q - —C G C F — 39 C F H— —C F T F n q 

y o o 


121 


9 


9 

+ • • • j 


C% - C G Tpn q + 66 C G Cp + 72 C\ - 24 C F T F n q 


(4.28) 


(4.29) 


(4.30) 


(4.31) 


where the ellipses refer to contributions that are lower order in g 3 . From eqs (4.10) and (4.28), I 
find the leading QCD 4-loop contribution to /3\: 


d 4) 


= yfatNcCp 


/ 3760 

w 


243094 


+C G T F n q 32^3 — 


( ^Cs - 1040Cs - 

5288 128tt 4 \ „ „ /7168 


243 


+ 


45 


243 
+ C G C F 


!|1) +Cc T f (H^ + 704 C 3 - 320( 5 


6608 , , 32t r 4 

— 6 + 14406 + is - 


o / _ 5884 327r 4 

+C F ( 96^3 + 320^5- 77 -1-— 


512 


-^ 2„2 


-T F n q — + T F n q ^ 


/21824 256 


3 ^ 


100 ^ ^ ^ ^ f 1124 , 64 32vr 4 

— 128 C F T F + C F T F n q I — f- — C 3 -g— 


+ 


(4.32) 


Now taking the limit e —>• 0, the effective potential is obtained in a loop expansion as 


y eff = 5^ 


£=0 


(167T 


2 \i 


v®. 


(4.33) 
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Note that unlike the loop expansion with bare parameters, eq. (3.1), here loop factors 1 /( 167 t 2 )^ 
have been extracted, similarly to eqs. (4.7) and (4.8). In terms of t and ln(i) defined in eqs. (1.2) 
and (1.3), the previously known results for the leading QCD effective potential contributions are: 


y(°) 

yd) 

y( 2 ) 


m 2 2 A 4 

T*? + 4 ^’ 

-N c t 2 [hi(f) — 3/2]. 


= g 2 N c C F t 2 [bln ( t ) — 161n(t) + 18]. 


(4.34) 

(4.35) 

(4.36) 


from ref. [22], and the three-loop result [24]: 


yd) = glN c C F t 2 {c G 

2609 44 4 232 


22— 3 , , 185— 2/ , , _ 1111,=-. , 

- xyIn (t) + —In (t) + (24( 3 -x—)hi(t) 

3 3 6 


+ -jy + y ^ 4 - + y ln 2 ( 2 )[vr 2 - ln 2 ( 2 )] - 128a 4 


+C F 


- 241n 3 (t) + 631n 2 (i) - (48C 3 + ^)ln(t) + y - y vr 4 


32 


+192C 3 - ^ ln 2 ( 2 )[ 7 r 2 - In 2 ( 2 )] + 256a 4 

O 


+T F n q 


+T f 


_ 232 

481n(f)- — + 96C 3 


(4.37) 


The new 4-loop result (with group-theory quantities left general) takes the form: 


=glC F N c t 2 Y J Y, <?ln n (i)y (4) (S,n), 

Q n =0 


(4.38) 


in terms of the group theory invariants in the set Q from eq. (3.13). The list of 50 coeffi¬ 
cients y) 4 )(l/,n) is again rather lengthy, and so is provided in another ancillary electronic file 
V4MSbar.txt. After substituting in the Standard Model values for the group theory constants, the 
result combines and simplifies to: 


yd) = 


git 2 


13820381 1747112C 3 1984C 5 40288(7 


270 


5888 

8 T 

+ln (t) 


ir 2 In 3 ( 2 ) - 

27680C 3 


45 

36064 


9 


298894vr 4 1780vr 6 5888 In 5 (2) 


1215 


405 
63200C 5 
9 


7T 4 In(2) + 


78464 


81 

1547146 

27 


In 2 (2) [In 2 (2) — 7 r 2 ] + 


243 

627712o 4 

27 


135 
47104o 5 
9 


2087T 4 640 9 

—— + — hr 2 (2) [In 2 (2) - vr 2 ] + 5120o 4 
y o 


+(30584 - 2400C 3 )ln 2 (f) - 91441n 3 (t) + 13801n 4 (f) 


(4.39) 


Equation (4.39) can be consistently added to the 3-loop effective potential as given in refs. [22] and 
[24], Also, the condition for the minimum v = r/> m j n of the Landau gauge effective potential of the 
Standard Model (including the effects of resummation of the Goldstone boson contributions from 
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the terms up to 3-loop order) is obtained by subtracting 


(167T 


2^4 


a 4 = 


i i dv ( 4 ) 


(167T 2 ) 4 V d(j) 


(4.40) 


computed using eq. (1.2) and (1.3) above, from the right-hand side of eq. (4.18) of ref. [25]. 


V. DISCUSSION 

The main results of this paper are the leading QCD 4-loop contributions to the Higgs self¬ 
coupling beta function f5\ and to the effective potential and its minimization condition. In each case, 
it is certainly possible that other contributions at 4-loop order, and the presently unknown 3-loop 
effects involving electroweak couplings in the case of the effective potential, could be numerically 
comparable to or even larger than the ones found here. The same is certainly true of parametric 
uncertainties from the top-quark Yukawa coupling (or mass) and the strong coupling. Therefore 
the results found here are perhaps most useful, for the present, as ways of formalizing estimates of 
purely theoretical error. 

The 4-loop leading QCD contribution of eq. (4.32) to the A beta function can be expressed in 
numerical form as 


A/ 3a = 8308.17^3^. (5.1) 

This can be compared to the leading QCD 1, 2, and 3-loop contributions: 

A '" 1 " 8 QCD = ^2 (-erf) + p^(- 329 ?rf) + pT_(— 100.40293V). ( 5 - 2 ) 

We see that the 4-loop contribution has a sign opposite to that of the other terms, and is larger 
in magnitude than one might have expected from a simple geometric progression. However, its 
magnitude is still only half as big as the 3-loop term in eq. (5.2) even at Q = M t , and in absolute 
terms it makes only a tiny difference in extrapolating A to high energy scales. 

The effective potential contribution of eq. (4.39) can similarly be expressed in numerical form 
as: 


U (4) = g®t 2 [ 5 93 66 .97 — 54056.36 ln(t) + 27699.06 ln 2 (t) — 91441n 3 (t) + 13801n 4 (i)l. (5.3) 


It follows that the corresponding contribution to the effective potential minimization condition 

(5.4) 


m + '-" r 2\2^ 2 (q 67r 2^ 3 n^2U^ 4 + ‘ ‘ ‘ 


167T 2 


(167T 2 


(167T 2 
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is, numerically: 

A 4 = g\y 1 2 t t [64677.58 - 52714.59 hi(t) + 27966.13 Ii) 2 (t) - 12768IK 3 (t) + 2760h) 4 5 6 (f)l , (5.5) 

where for £ = 1, 2,3 were given in ref. [25]. Consider the VEV and other MS parameters of the 
Standard Model at benchmark values 


v(M t ) = 246.647 GeV, (5.6) 

A (M t ) = 0.12597, (5.7) 

y t (M t ) = 0.93690, (5.8) 

g 3 (M t ) = 1.1666, (5.9) 

g(M t ) = 0.647550, (5.10) 

g \M t ) = 0.358521, (5.11) 


at Q = M t = 173.34 GeV. These choices provide agreement with the measured values of the h, W, 
and Z boson masses in the pure MS scheme [67-69]. Using only the previously known 3-loop con¬ 
tributions in eq. (5.4), the resulting Higgs squared mass parameter is: m 2 (M t ) = —(92.890 GeV) 2 . 
Now including the new contribution of eq. (5.5) gives instead m 2 (M t ) = —(92.926 GeV) 2 . Thus I 
find 

A (V- m 2 ) = 36 MeV (5.12) 

from the leading QCD 4-loop contribution, at the scale Q = Mf. The parameter m 2 is not directly 
constrained by experiment, but it can be connected to ultraviolet completions that may predict 
it in terms of other underlying parameters that can be measured, eventually. This could occur in 
models of supersymmetry breaking, for example. 
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